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ABSTRACT 


Procedures  are  presented  whereby  natural  vibration 
and  transient  response  analyses  of  partially  submerged, 
supercavitating  marine  propellers  may  be  calculated  by 
means  of  a  standard  NASTRAN  structural  analysis  program. 

The  use  of  these  procedures  was  demonstrated  by  analyses 
of  an  existing  propeller,  for  which  some  experimental  data 
were  available.  Results  were  comparable  only  in  a  limited 
way,  but  peak  calculated  dynamic  stresses  were  of  the  same 
order  of  magnitude  as  those  observed  experimentally,  and 
were  apparently  more  than  twice  the  magnitude  of  equivalent 
static  stresses . 


ADMINISTRATIVE  INFORMATION 

Funds  for  this  investigation  were  obtained  from  the  Naval  Material 
Command  (NAVMAT  0333)  in  support  of  an  ongoing  Ship  Performance  and  Hydro¬ 
mechanics  Exploratory  Development  Program  (Program  Element  62543N,  Task 
Area  ZP-43-421-001)  assigned  to  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC).  The  High  Speed  Propulsor  Task  of  this 
Exploratory  Development  Program  (DTNSRDC  Work  Unit  No.  1300-200)  provided 
funds  to  complete  work  presented  in  this  report. 

INTRODUCTION 

The  utilization  of  high  speed  marine  vehicles  by  the  United  States 
Navy  has  created  a  need  for  an  improved  structural  analysis  and  design 
capability  for  high  speed  propulsors  .  These  are  generally  supercavitating 
propellers,  which  may  also  be  partially  submerged.  In  their  hydrodynamic 
performance  and  structural  integrity,  these  propellers,  because  of  their 
different  geometries  and  modes  of  operation,  represent  substantial  depar¬ 
ture  from  conventional  propellers,  for  which  a  large  body  of  experience 
exists.  Dynamic  aspects  of  structural  response,  in  particular,  are  ex¬ 
pected  to  be  more  significant  tor  these  propellers.  These  aspects  need  to 
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be  evaluated  so  that  attempts  to  attain  optimum  hydrodynamic  designs  will 
not  be  limited  unnecessarily  by  structural  considerations.  An  effort  was 
therefore  instituted  to  develop  methods  of  predicting  dynamic  stress  and 
deflection  behavior  of  high  speed  propellers  under  realistic  time-dependent 
loading  conditions. 

A  preliminary  investigation  indicated  that  meaningful  results  were 
not  likely  to  be  obtained  other  than  by  computer  studies.  Various  computer 
programs  have  been  used  to  calculate  the  full  vibration  modes  and  frequen¬ 
cies  of  completely  submerged  propellers,  but  until  recently  the  effect  of 
immersion  in  water  has  been  accounted  for  by  assuming  that  a  more  or  less 
arbitrary  mass  of  fluid  is  attached  to  the  propeller  blade.  This  "added 
effective  mass"  approach,  which  incorporates  the  assumption  that  the  fluid 
is  incompressible,  is  considered  to  account  adequately  for  fluid  structure 
interaction  in  most  vibration  and  transient  response  problems,  provided 
valid  added  fluid  mass  terms  can  be  obtained.  A  few  investigators  have 
recently  published  descriptions  of  analyses  of  dynamic  behavior  of  sub¬ 
merged  structures  in  which  effects  of  submergence  in  water  are  computed  on 
a  rational  basis.  One  such  investigation^*  yielded  a  procedure  by  which 
the  existing  NASTRAN  finite  element  structural  analysis  program  could  be 
applied  to  the  submerged  structure  problem,  even  without  making  the  incom¬ 
pressibility  assumption.  That  work  forms  the  basis  of  the  analyses  de¬ 
scribed  here.  The  approach  was  selected  because  of  the  availability  of 
NASTRAN  at  DTNSRDC  and  because  of  its  flexibility  of  application. 


*A  complete  listing  of  references  is  given  on  page  57. 
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Further  modifications  were  made  to  this  procedure  to  adapt  it  to  the 
general  problem  of  predicting  free  vibration  frequencies  and  transient 
response  of  supercavitating ,  partially  submerged  propellers.  The  resulting 
approach  was  applied  to  determine  natural  frequencies  of  two  propellers, 
and  the  dynamic  stress  response  of  one  propeller  for  which  some  experimen¬ 
tal  data  were  available.  Only  a  limited  comparison  of  results  was  possible; 
however,  the  calculated  frequencies  and  dynamic  stresses  appeared  to  be 
reasonable.  The  dynamic  analysis  and  its  application  via  NASTRAN  are 
described  in  this  report. 


DESCRIPTION  OF  ANALYSIS 

In  the  analysis  to  be  used  here,  the  fluid  surrounding  the  submerged 
structure,  in  this  case  an  individual  propeller  blade,  is  represented  by  a 
finite  element  idealization.  More  precisely,  a  quantity  of  fluid  suffi¬ 
ciently  large  that  its  effect  on  the  structure  should  nearly  be  the  same 
as  that  of  an  infinite  quantity  is  represented.  The  theory  governing  this 
approach  is  described  by  Zienkiewicz . ^  A  concise  presentation,  with  a 
delineation  of  the  strategy  of  application  by  means  of  existing  NASTKAN 
capabilities  is  given  by  Everstine  et  al.^  The  theory  and  application 
procedure  will  be  briefly  summarized  here. 

Conventional  analytical  techniques  lead  to  the  following  matrix  equa¬ 
tion,  which  describes  the  dynamic  behavior  of  the  propeller  blade  when 
represented  in  finite  element  form: 

[Ml  {u  }  +  [  K  ]  { 11 1  =  IF}  +  [L](pf  (1) 
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where  { u }  is  the  set  of  unknown  grid  point  displacements  and  {pi  is  a  set 
of  applied  time-dependent  forces.  The  column  matrix  Cp}  is  the  set  of 
unknown  pressures  on  the  fluid-structure  interface  grid  points.  The  last 
term,  [ L ] { p }  »  represents  the  fluid  loading  on  the  structure. 

The  set  of  unknown  pressures  at  all  the  grid  points  of  the  fluid  fi¬ 
nite  element  idealization  is  represented  in  the  following  matrix  equation: 

[Q]  ip  ■  +  [H]  ip  'f  =  -  [S]  iu  •  (2) 

The  last  term  in  Equation  (2)  represents  the  effect  of  accelera t ions  of  the 
structural  grid  points  at  the  fluid  structure  interface  on  the  pressures 
in  the  fluid.  The  bar  placed  over  iu  •  indicates  that  only  normal  interface 
deflections  are  involved  in  that  term.  The  fluid  may  be  considered  an 
acoustic  medium,  in  which  case  Equation  (2)  provides  a  valid  description 
of  its  dynamic  behavior  if  it  actually  represents  the  wave  equation 


with  appropriate  boundary  conditions.  In  Equation  (3),  c  represents  the 
speed  of  sound  in  the  fluid.  It  was  shown  previously*  that  the  implemen¬ 
tation  of  the  procedure  can  be  accomplished  by  means  of  a  standard  struc¬ 
tural  program.  The  left  side  of  Equation  (2)  is  treated  as  being  analogous 
to  the  left  side  of  Equation  (1),  i.e.,  it  is  treated  as  a  dynamics  problem 
with  unknown  variables  tpl  .  Only  one  translational  degree  of  freedom  at 
each  grid  point  of  the  fluid  is  left  unconstrained,  however,  and  that  one 
remaining  component  is  interpreted  as  pressure.  In  addition,  all  degrees 
of  freedom  may  be  constrained  on  the  outer  fluid  boundary  (all  boundary 
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points  not  in  contact  with  the  moveable  structure).  Finally  for  the 
three-dimensional  case,  the  following  material  property  matrix,  repre¬ 
senting  Hooke's  law,  is  used: 
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With  the  NASTRAN  program,  it  is  not  usually  necessary  to  input  the  matrix 
(4).  As  pointed  out  by  Everstine,  the  same  information  can  be  input  by 
specifying  an  isotropic  material  with  the  following  values: 

Shear  modulus  <J  =  1.0 

Poissons  ratio  v  is  not  input  (so  the  program  will 
compute  a  value). 

Young's  modulus  E  =  *0;  x  '■'*  1.0 

Here,  a  +  l.u  should  be  numerically  indistinguishable  from  >.  On  most 
computers,  a  =  10^  will  suffice.  Equations  (1)  and  (2)  may  now  be  com¬ 
bined  into  one  set: 


M 

0 

5 

♦ 

K 

-L 

r 

= 

F 

s 

Q. 

p 

0 

H  _ 

Ip 

0 

(5) 


A  dynamic  problem  governed  by  Equation  (.5)  can  be  run  (the  free  vibra¬ 
tion  mode  analysis  results  from  letting  F  be  null)  by  providing  finite  ele¬ 
ment  data  for  both  structure  and  fluid,  in  which  case  all  matrices  except 


S  and  L  in  Equatio 


(5)  will  be  set  up  automatically.  Note  that  grid 


points  for  structure  and  fluid,  including  those  on  the  interface  which  have 
the  same  coordinates,  are  numbered  separately.  Thus,  the  only  coupling  be¬ 
tween  fluid  and  structure  is  provided  by  the  off-diagonal  matrices,  S  and  L. 
Equation  (1),  which  expresses  a  balance  of  forces,  shows  that  L  relates  the 
pressure  at  the  fluid-structure  interface  to  force  resulting  from  that 
pressure.  In  a  lumped  formulation,  the  terms  of  L  may  be  taken  to  be  the 
interface  surface  areas  associated  with  each  fluid-structure  interface  grid 
point.  In  general,  there  will  be  three  areas  for  each  grid  point  -  one 
for  each  coordinate  direction.  Zienkiewic2^  has  shown  that 

S  =  >'L F  (6) 

Equation  (5)  then  becomes 


T 

The  L  and  L  matrices  may  be  determined  separately  and  input  directly, 
using  the  UMIG  option  of  NASTRAN.  Thus,  formulation  of  Equation  (7)  (or 
Equation  (5))  is  completed,  and  the  solution  is  carried  out  by  the  program. 
However,  the  presence  of  the  off-diagonal  matrices  destroys  the  symmetry 
of  the  mass  and  stiffness  matrices,  requiring  the  use  of  time-consuming 
complex  eigenvalue  routines  in  the  solution  of  these  equations. 

The  natural  vibration  frequencies  for  some  flat  plates  with  varying 
levels  of  immersion  were  calculated  by  the  method  described  and  compared 
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with  experimental  results  by  Marcus.^  Good  agreement  was  obtained  when  a 
sufficient  amount  of  water  was  modeled,  and  in  fact  the  results  provide 
a  useful  reference  for  judging  the  size  of  a  fluid  mass  required  to  repre¬ 
sent  an  infinite  body  of  water. 

APPLICATION  TO  PROPELLER  BLADES 

It  was  planned  originally  to  apply  the  analysis  as  represented  by 
Equation  (7)  to  a  supercavitating  propeller  blade,  P3604,  for  which  data 
for  a  finite  element  idealization  were  available  from  a  previous  stress 
study.  The  idealization  consisted  of  155  grid  points  and  16  quadratic 
isoparametric,  twenty-noded ,  three-dimensional  elements  (NASTRAN  data  card 
CIHEX2).  The  purpose  of  this  sample  calculation  was  to  verify  that  appli¬ 
cation  of  this  procedure  to  the  more  complicated  propeller  blade  would  be 
feasible.  It  was  apparent,  however,  that  the  cost  of  such  a  computer  run 
could  be  high.  A  fluid  idealization  was  constructed  of  only  three  layers 
of  elements  and  extended  one  element  beyond  the  edges  of  the  blade.  The 
fluid  idealization,  which  was  thought  to  be  just  large  enough  to  represent 
an  infinite  quantity,  but  only  on  one  side  of  the  blade,  contained  90 
three-dimensional  elements  and  578  grid  points.  A  total  of  263  uncon¬ 
strained  degrees  of  freedom  resulted,  and  this  number  would  be  about 
doubled  for  fluid  on  both  sides  of  the  blade.  The  blade  itself  was  con¬ 
sidered  fixed  along  its  base,  so  that  a  total  of  396  degrees  of  freedom 
(3  for  each  unconstrained  grid  point)  remained.  The  total  number  of 
degrees  of  freedom  for  fluid  and  structure  thus  led  to  mass  and  stiffness 
matrices  for  Equation  (7)  of  order  659.  Eigenvalue  extraction  of  these 
large  unsymmetrical  matrices  would  have  required  a  time-consuming  complex 


extraction  routine.  It  was  decided,  therefore,  to  attempt  to  reduce  the 


computer  effort  required  without  any  significant  sacrifice  of  accuracy. 
In  Equations  (2)  and  (7),  the  Q  matrices,  which  contain  a  factor  of 


1/c,  represent  the  effect  of  compressibility  of  the  fluid.  This  effect  is 
not  considered  significant  except  for  high  frequencies  or  in  confined 
fluid  spaces.  Most  analyses,  including  that  of  Everstine,*  have  incor¬ 
porated  the  simplifying  assumption  of  incompressibility;  as  a  result,  1/c, 
and  hence  Q,  are  zero.  This  assumption  was  not  included  in  the  formulation 
of  Equation  (7),  which  was  set  up  for  ready  computer  application,  because 
there  was  little  advantage  in  doing  so. 

When  Q  is  taken  to  be  zero,  Equation  (7)  is  readily  reduced  to 

[M+f>LH~1LT]  fu  }  +  [ K ]  f u  }  =  { F  ;  (8) 

One  additional  significant  operation,  that  of  obtaining  the  inverse 

of  H,  is  required  to  set  up  Equation  (8).  However,  this  and  the  eigenvalue 

extraction  operations  are  performed  on  matrices  about  one-half  the  size 

of  those  of  Equation  (7).  Furthermore,  the  matrices  of  Equation  (8)  are 
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symmetrical.  The  term  pLH  L  is  interpreted  as  the  added  effective  mass 
of  the  fluid.  It  is  clear  that  fluid  on  only  one  side  of  the  propeller 
blade  may  be  idealized,  and  the  added  mass  term  can  be  made  to  represent 
fluid  on  both  surfaces  by  using  a  value  for  fluid  density,  p,  that  is 
double  the  actual  value. 

No  standard  NASTRAN  solution  option,  or  "rigid  format,"  is  available 
to  formulate  and  solve  Equation  (8).  Executive  modules  can  be  manipulated, 
however,  though  the  Direct  Matrix  Abstraction  Program  (DMAP)  capability, 
and  existing  rigid  formats  may  be  altered  in  t lie  same  way.  These  capa¬ 
bilities  were  used  to  devise  a  procedure  for  implementing  Equation  (8). 


First,  the  fluid  matrix  (or  rather,  (l/p)H)  was  computed  and  made  available 
for  later  reference.  This  fluid  "stiffness"  matrix  was  produced  by  making 
a  static  stress  run  on  the  fluid  mass  with  material  properties  and  boundary 
conditions  as  discussed  above.  ASET1  cards  were  used  to  specify  that  the 
solution  set  be  retained  only  at  the  interface  grid  points,  and  an  "0UTPUT1" 
DMAP  ALTER  instruction  to  specify  that  the  resulting  condensed  stiffness 
matrix  (or  H)  be  saved  on  a  permanent  file.  Another  ALTER  instruction  was 
used  to  stop  further  execution  of  the  problem.  This  run  can  also  be  used 
to  obtain  terms  of  the  L  matrix.  It  was  suggested  previously  that  areas 
associated  with  each  structural  interface  grid  point  can  be  used  as  L 
matrix  terms .  Calculating  these  areas  and  obtaining  their  projections  in 
the  three  coordinate  directions  for  a  propeller  surface  can  be  quite 
difficult.  It  is  much  easier,  and  equivalent,  to  specify  a  unit  pressure 
loading  on  the  interface  and  request  a  printout  of  the  three  equivalent 
nodal  forces  on  these  grid  points.  The  forces  will  thus  be  calculated  in 
a  consistent  manner.  These  forces  can  be  used  to  make  up  the  L  matrix. 

They  are  functions  of  only  the  pressurized  surface  shape  and  the  location 
of  the  grid  points.  Therefore,  since  these  two  quantities  are  identical 
for  both  the  structure  and  the  fluid  at  their  interface,  the  unit  pressure 
load  can  be  applied  to  the  idealized  fluid  interface  in  the  same  run  from 
which  the  H  matrix  is  obtained.  The  equivalent  nodal  forces  can  then  be 
entered  on  DMIG  cards  as  L  matrix  terms.  Since  the  fluid  run  has  only  one 
degree  of  freedom,  while  all  three  translational  degrees  are  needed  to 
obtain  the  three  components  for  each  surface  grid  point,  the  equivalent 
force  matrix  must  be  printed  out  before  constraints  are  applied.  This 
printout  is  requested  through  another  ALTER  instruction. 
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To  summarize  the  procedure:  a  static  stress  run  is  started  on  the 
fluid  mass  with  unit  pressure  applied  to  the  interface.  All  interface  grid 
points  are  specified  on  ASET1  cards.  The  stiffness  matrix  (after  con¬ 
straints  and  condensation  are  applied)  is  placed  on  a  permanent  file  via 
the  required  control  card.  An  0UTPUT1  ALTER  instruction  and  the  equivalent 
nodal  forces  are  then  printed  out,  at  which  point  the  calculation  is 
stopped.  The  matrix  condensation  is  not  theoretically  required,  but  is 
done  to  save  space;  it  does  not  affect  the  final  results. 

Equation  (8)  may  be  set  up  for  the  propeller  and  solved  to  obtain  ei¬ 
ther  the  frequency  response  or  the  dynamic  response  to  applied  time-depen¬ 
dent  forces  by  running  the  propeller  data  with  the  additional  data  already  • 
obtained.  The  natural  frequencies  are  obtained  by  specifying  a  modified 
natural  modes  calculation  with  F  null  (Rigid  Format  3);  the  dynamic 
response  is  obtained  by  specifying  a  transient  response  calculation  with 
F  input  (Rigid  Format  9).  Each  of  these  rigid  formats  requires  its  own 
ALTER  routine,  although  each  routine  works  in  essentially  the  same  way, 
and  the  additional  input  is  handled  similarly.  Structural  grid  point  data 
must  include  one  extra  scalar  point  (SPOINT)  for  each  fluid  interface  grid 
point  listed  on  ASET1  cards  and  thus  retained  as  a  degree  of  freedom  in 
the  H  matrix  calculation.  These  scalar  points  are  dummy  grid  points  which 
have  no  coordinates,  but  have  one  allowable  degree  of  freedom  each  (the 
zero  component).  They  are  established  by  listing  their  identification 
numbers  on  SPOINT  cards.  These  identification  numbers  are  (for  conven¬ 
ience)  the  fluid  interface  grid  point  numbers,  and  they  must  be  different 
from  all  propeller  grid  point  numbers.  The  purpose  of  the  scalar  points 
is  to  expand  the  dimensions  of  t  tie  structural  stiffness  matrix,  belore 
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constraints  are  applied,  sufficiently  to  accept  the  L  (or  L  )  matrix  which 
is  also  input  to  the  stiffness  matrix  via  DM1G  cards.  It  was  found  more 
advantageous  to  require  the  L  matrix;  the  L  matrix  can  readily  be  derived 
from  it  by  transposition.  The  terms  of  the  matrix  are  specified  on  DM1G 
cards,  each  of  which  relates  one  of  the  three  nodal  force  components  per 
unit  pressure  at  a  wetted  propeller  grid  point,  obtained  from  the  fluid 

calculation  with  the  zero  component  of  the  corresonding  scalar  point.  If 

T 

the  algebraic  signs  of  all  of  the  L  terms  are  reversed,  those  of  L  will 

—  1  T 

also  be  reversed,  and  the  added  mass  term,  pLH  L  ,  will  be  unchanged; 

T 

therefore  a  sign  error  for  the  L  matrix  as  a  whole  has  no  effect.  Since 
only  the  interface  fluid  points  are  input  here  (as  scalar  points),  the  zero 
rows  of  L  are  eliminated  so  that  the  matrix  conforms  with  the  condensed 
version  of  H  .  The  (l/p)H  matrix  is  obtained  from  the  file  on  which  it  has 
been  stored  by  including  the  required  control  card  and  an  INPUTT1  ALTER 
instruction. 

Two  partitioning  matrices,  named  PVEC1  and  PVECD1 ,  must  also  be  input 

via  DMI  cards.  These  are.  column  matrices  having  the  dimensions  (N+S)xl.  N 

is  the  total  number  of  structural  degrees  of  freedom  (number  of  grid  points 

x6)  and  S  is  the  number  of  scalar  points.  For  PVEC1 ,  the  first  N  terms 

are  zero,  and  the  remaining  S  terms  are  1.0.  PVEC1  is  used  to  partition 

T 

the  stiffness  matrix  to  obtain  L  ,  used  with  (l/p)H  by  the  ALTER  routine  to 

-1  T 

form  the  added  mass  term  pLH  L  ,  which  is  then  added  to  the  structural 
mass  matrix.  Note  that  the  inverse  of  (l/p)H  is  pH  ^ .  For  the  PVEDC1  matrix, 
all  (N+S)  terms  are  zero.  PVECD1  is  used  to  form  a  "dummy"  partition  on 
tiie  "augmented"  mass  matrix,  which  now  includes  the  added  mass  terms,  so 
that  tlie  program  will  correctly  perceive  this  matrix  as  symmetrical. 
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One  final  type  of  additional  input  is  required:  all  scalar  points  must 
be  listed  on  single  point  constraint  (SPC)  cards.  These  constraints  are  in 
addition  to  whatever  constraints  are  applied  to  the  propeller.  It  is  sug¬ 
gested  that  20-node  three-dimensional  CIHEX2  elements,  with  two  integration 
points  specified,  be  used  to  represent  both  the  fluid  and  the  propeller, 
and  that  ttie  consistent  mass  matrix  formulation  be  specified  for  the  pro¬ 
peller. 

The  procedure  just  described  can  be  used  to  obtain  vibration  modes 
and  dynamic  response  of  propellers  under  a  variety  of  conditions  of  sub¬ 
mergence.  Partial  submergence  at  any  given  instant  can  be  represented 
during  a  fluid  run  by  listing  on  ASET1  cards  only  the  fluid  grid  points 
which  correspond  to  submerged  points  on  the  propeller.  During  a  structural 
run  these  points  are  listed  on  SPOINT  cards.  Also,  the  dimensions  of  the 
PVECi  and  PVECD1  column  matrices  must  be  adjusted  by  making  the  number  of 
terms  with  value  of  1.0  at  the  end  of  the  PVEC1  column  equal  to  the 
number  ot  wetted  grid  points.  Wetting  of  both  surfaces  of  the  blade  is 
simulated  by  using  double  the  actual  value  of  fluid  density.  •  ,  for  the 
fluid  run.  Wetting  of  only  one  surface,  such  as  occurs  during  supercavi¬ 
tation,  is  represented  by  use  of  the  actual  value  of  ,  . 

This  procedure  can  be  extended,  without  further  DMAP  program  modifi¬ 
cations,  to  analyze  dynamic  response  of  partially  submerged  propellers 
and  to  take  into  account  the  progressively  changing  degree  of  submergence. 
The  extended  procedure  consists  essentially  of  running  the  analysis  in 
separate  steps,  one  for  each  arbitrarily  chosen  stage  of  submergence.  A 
convenient  way  to  determine  a  new  stage  oi  submergence  is  to  consider  the 
wetting  ot  a  new  row  of  propeller  elements.  This  extended  procedure  was 
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devised  and  was  applied  to  an  existing  propeller  for  which  some  experi¬ 
mental  and  analytical  data  were  available.  The  details  of  this  application 
will  be  presented  later.  The  finite  element  idealization  of  this  blade 
was  laid  out  so  that  one  series  of  element  boundaries  coincided  with  water 
lines  on  the  propeller  blades  as  it  progressed  into  the  water  during  each 
revolution.  There  were  two  such  boundaries  in  this  rather  coarse  ideali¬ 
zation,  defining  three  rows  of  elements;  hence,  three  stages  of  submergence 
were  assumed. 

The  fluid  run,  to  obtain  the  condensed  (l/p)H  matrix,  was  run  once  for 

each  stage  of  submergence.  For  each  run,  the  solution  set  of  wetted  grid 

points,  as  defined  by  ASET1  cards,  contained  only  the  points  presently 

wetted,  i.e.,  points  just  submerged  plus  all  points  previously  submerged. 

Similarly,  unit  pressure  was  applied  to  the  corresponding  wetted  element 

T 

faces  to  obtain  the  L  matrix  terms  for  the  same  wetted  points.  Finally, 
those  interface  grid  points  not  yet  submerged  had  to  be  considered  part 
of  the  fluid  boundary  and  were  constrained.  Constraints  were,  of  course, 
removed  from  those  points  as  they  became  submerged  in  subsequent  stages  and 
added  to  those  included  on  ASETl  cards.  Each  of  the  (l/p)H  matrices  thus 
produced  and  stored  on  a  permanent  file  was  assigned  a  separate  name;  in 
this  case  the  names  used  were  HMATRIXSTEP 1 ,  HMATR1XSTEP2 ,  and  HMATRIXSTEP3 . 
Such  separate  identification  of  the  L  matrices  was  not  necessary,  as 
they  are  manually  punched  on  UMIG  cards  and  input  separately  to  the  struc¬ 
tural  computer  runs.  The  structural  run  in  this  case  was  also  carried 
out  once  for  each  stage  of  submergence,  but  all  except  the  first  were 
restart  runs.  Thus,  all  but  the  last  run  was  checkpointed. 
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The  series  of  runs  is  set  up  as  follows.  First,  time-dependent 
forces  are  specified  for  the  selected  grid  points  by  means  of  the 
required  DAREA,  DLOAD,  TABLED1 ,  and  TLOADi  cards.  The  time  spanned  by 
these  forces  covers  the  entire  time  for  which  the  analysis  is  applied. 

The  time  at  which  the  force  is  initiated  at  each  selected  grid  point  is 
dictated  by  the  speed  of  rotation  of  the  propeller  and  by  the  location  of 
the  point.  The  user  must  also  specify  the  length  of  each  integration  time 
step,  and  the  number  of  steps,  on  a  TSTEP  card.  The  calculation  is  ter¬ 
minated  when  the  time  defined  by  the  product  of  time  step  x  number  of  steps 
on  the  TSTEP  card  has  elapsed.  If  a  uniform  time  step  is  used  throughout, 
the  number  of  steps  may  be  specified  to  make  each  run  cover  just  the  length 
of  time  required  for  the  water  line  to  progress  from  the  first  to  the 
second  boundary  defining  the  present  stage  of  submergence.  The  first  run 
and  each  restart  then  require  a  different  TSTEP  card  (unless,  by  coinci¬ 
dence,  each  stage  of  submergence  lasts  the  same  length  of  time).  Other 
changes  must  also  be  made  between  runs.  Scalar  points  for  all  interface 
fluid  grid  points  which  are  to  be  wetted  at  any  time  during  the  analysis 
must  be  input  (on  SP01NT  cards)  for  the  first  run,  and  this  does  not  change 
for  succeeding  restarts.  The  same  is  true  of  the  scalar  point  constraints. 
The  PVECDl  matrix  is  also  input  as  before,  with  column  dimension  equal  to 
six  times  the  number  of  propeller  grid  points  plus  the  number  of  scalar 
points,  and  it  is  not  changed.  The  PVEC1  column  matrix  is  now  replaced 
by  a  series  of  matrices,  one  for  each  structural  run,  and  they  may  be 
designated  PVEC1  ,  PVEC2 ,  and  PVEC3 .  The  column  dimension  is  the  same  for 
all  of  these,  and  is  the  same  as  that  of  PVECDl.  The  number  and  location 
ot  the  "l.U"  terms,  however,  must  correspond  to  the  total  set  of  wetted 
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devised  and  was  applied  to  an  existing  propeller  for  which  some  experi¬ 
mental  and  analytical  data  were  available.  The  details  of  this  application 
will  be  presented  later.  The  finite  element  idealization  of  this  blade 
was  laid  out  so  that  one  series  of  element  boundaries  coincided  with  water 
lines  on  the  propeller  blades  as  it  progressed  into  the  water  during  each 
revolution.  There  were  two  such  boundaries  in  this  rather  coarse  ideali¬ 
zation,  defining  three  rows  of  elements;  hence,  three  stages  of  submergence 
were  assumed. 

The  fluid  run,  to  obtain  the  condensed  (l/p)H  matrix,  was  run  once  for 
each  stage  of  submergence.  For  each  run,  the  solution  set  of  wetted  grid 
points,  as  defined  by  ASET1  cards,  contained  only  the  points  presently 
wetted,  i.e.,  points  just  submerged  plus  all  points  previously  submerged. 
Similarly,  unit  pressure  was  applied  to  the  corresponding  wetted  element 
faces  to  obtain  the  L  matrix  terms  for  the  same  wetted  points.  Finally, 
those  interface  grid  points  not  yet  submerged  had  to  be  considered  part 
of  the  fluid  boundary  and  were  constrained.  Constraints  were,  of  course, 
removed  from  those  points  as  they  became  submerged  in  subsequent  stages  and 
added  to  those  included  on  ASET1  cards.  Each  of  the  (1/ (>) H  matrices  thus 
produced  and  stored  on  a  permanent  file  was  assigned  a  separate  name;  in 
this  case  the  names  used  were  HMATR1XSTEP1 ,  HMATK1XSTEF2  ,  and  HMATR1XSTEP  . 

t 

Such  separate  identification  of  the  L  matrices  was  not  necessary,  as 
they  are  manually  punched  on  DM1G  cards  and  input  separately  to  the  struc¬ 
tural  computer  runs.  The  structural  run  in  this  case  was  also  carried 
out  once  for  each  stage  of  submergence,  but  all  except  the  first  were 
restart  runs.  Thus,  ail  but  the  last  run  was  checkpointed. 


The  series  of  runs  is  set  up  as  follows.  First,  time-dependent 
forces  are  specified  for  the  selected  grid  points  by  means  of  the 
required  DAREA,  DLOAD,  TABLED1 ,  and  TLOADi  cards.  The  time  spanned  by 
these  forces  covers  the  entire  time  for  which  the  analysis  is  applied. 

The  time  at  which  the  force  is  initiated  at  each  selected  grid  point  is 
dictated  by  the  speed  of  rotation  of  the  propeller  and  by  the  location  of 
the  point.  The  user  must  also  specify  the  length  of  each  integration  time 
step,  and  the  number  of  steps,  on  a  TSTEP  card.  The  calculation  is  ter¬ 
minated  when  the  time  defined  by  the  product  of  time  step  x  number  of  steps 
on  the  TSTEP  card  has  elapsed.  If  a  uniform  time  step  is  used  throughout, 
the  number  of  steps  may  be  specified  to  make  each  run  cover  just  the  length 
of  time  required  for  the  water  line  to  progress  from  the  first  to  the 
second  boundary  defining  the  present  stage  of  submergence.  The  first  run 
and  each  restart  then  require  a  different  TSTEP  card  (unless,  by  coinci¬ 
dence,  eacii  stage  of  submergence  lasts  the  same  length  of  time).  Other 
changes  must  also  be  made  between  runs.  Scalar  points  for  all  interface 
tiuid  grid  points  which  are  to  be  wetted  at  any  time  during  the  analysis 
must  be  input  (on  SP01NT  cards)  for  the  first  run,  and  this  does  not  change 
for  succeeding  restarts.  The  same  is  true  of  the  scalar  point  constraints. 
The  PVECD1  matrix  is  also  input  as  before,  with  column  dimension  equal  to 
six  times  the  number  of  propeller  grid  points  plus  the  number  of  scalar 
points,  and  it  is  not  changed.  The  PVEC1  column  matrix  is  now  replaced 
by  a  series  of  matrices,  one  for  each  structural  run,  and  they  may  be 
designated  PVEC1,  PVEC2 ,  and  PVEC3 .  The  column  dimension  is  the  same  for 
all  of  these,  and  is  the  same  as  that  of  PVECDl  .  The  number  and  location 
ot  the  "1.U"  terms,  however,  must  correspond  to  t  he  total  set  of  wetted 


interface  points  as  represented  by  the  scalar  points  during  each  run.  The 
DMAP  ALTER  instruction  which  utilized  PVEC  must  also  be  changed  from  run 
to  run  to  call  for  the  currently  available  PVECi  matrix.  Another  ALTER 
instruction,  as  well  as  a  control  card,  must  be  changed  lor  each  run  to 
obtain  the  required  (l/p)H  matrix  for  each  stage.  Finally,  each  run  must 
also  have  a  new  L^  matrix  input  via  DM 1C  cards.  Since  this  information 
is  part  of  the  bulk  data,  the  previous  LT  matrix  must  be  eliminated  during 
each  restart. 

Sample  input  for  a  vibration  problem  and  for  a  transient  response 
problem  are  presented  in  the  appendix. 

VIBRATION  ANALYSIS  OF  P3604  PROPELLER 

The  vibration  analysis  procedure  was  lirst  applied  to  a  supercavitating 
propeller,  P36U4,  for  which  finite  element  idealization  data  were  already 
available  from  a  previous  stress  study.  The  blade  was  divided  into  lb 
isoparametric,  20-noded ,  three-dimensional  elements,  with  133  grid  points 
(Figure  1).  A  quantity  of  water  three  elements  thick,  extending  one  ele¬ 
ment  beyond  the  edges  of  the  blade,  and  generally  following  t  lie  twist  ol 
the  blade,  was  assumed.  The  elements  were  roughly  cubical,  and  the  thick¬ 
ness  of  the  water  modeled  was  about  one-half  to  two-thirds  the  length  ol 
the  blade.  This  size  was  judged  to  be  sufficient  to  represent  an  infinite 
body  of  water,  based  on  a  study  of  the  results  presented  by  Marcus. ^  With 
the  procedure  used  here,  larger  representations  could  be  investigated  at 
fairly  moderate  cost.  The  double  value  of  water  density  was  input  to 
represent  water  on  both  sides  ot  the  blade.  This  fluid  idealization  con¬ 
sisted  ol  37b  grid  points  and  90  elements.  A  total  ot  2(>3  degrees  ot 
freedom  remained  alter  constraints  were  applied.  The  fluid-structure 
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interface  surface  contained  65  grid  points,  which  were  retained  in  the 


T 

H  matrix.  The  L  matrix  terms  numbered  65x3  or  195.  However,  nine  of 

the  interface  points  were  adjacent  to  the  base  grid  points  of  the  propeller, 

T 

which  were  later  to  be  constrained.  The  corresponding  9x3,  or  27  L  matrix 

terms  were  therefore  not  included  on  DMIG  cards.  Space  was  later  reserved 

for  them  in  the  propeller  matrices,  however,  as  65  scalar  points  were 
T 

input,  but  the  L  boundary  terms  which  were  not  input  were  implicitly  zero. 
It  was  not  necessary  to  leave  out  these  terms;  they  would  be  eliminated 
later  by  constraints.  With  this  procedure,  the  pressures  in  the  water  at 
the  nine  base  grid  points  were  not  constrained  to  zero,  even  though  deflec¬ 
tions  of  the  corresponding  propeller  grid  points  were  constrained.  The 
tive  iowest  naturai  frequencies  of  the  blade  in  water  were  calculated  and 
are  presented,  with  the  corresonding  calculated  in-air  frequencies,  in 
Table  1.  The  modal  deflection  data  indicated  that  the  first  mode  of 
vibration  consisted  of  nearly  pure  spanwise  bending,  while  the  other  modes 
contained  a  combination  of  torsion  and  chordwise  bending.  Although  no 
experimental  data  exist  to  compare  with  this  sample  calculation,  the 
results  appear  to  be  consistent  with  experimental  results  for  more  conven¬ 
tional  blades,  particularly  for  the  first  mode.  It  seems  rather  remarkable 
that  the  ratios  of  frequencies  in  water  to  those  in  air  are  fairly  constant, 
ranging  from  0.46  to  0.55,  for  all  modes.  This  was  not  expected  because 
of  the  large  differences  in  thickness  of  the  blade  and  the  variety  ol 
possible  vibration  modes.  This  phenomenon  may  not  be  characteristic  of 
other  blade  geometries. 

The  version  of  NASTKAN  used  for  these  analyses  was  Level  17.  Proce¬ 
dures  for  both  the  vibration  and  transient  analyses  are  listed  In  the 
Append ix. 


16 


TRANSIENT  RESPONSE  ANALYSIS  OF  SES-IOOB  PROPELLER 


After  this  work  was  underway,  some  full  scale  test  results,  including 
both  vibration  and  transient  response  data,  became  available^’6  for  the  SES- 
IOOB  controllable  pitch,  partially  submerged,  supercavitating  propeller 
(number  701b. 3).  Allison's  work'1  also  contains  finite  element  analyses  of 
the  vibration  modes  and  transient  response  of  this  blade,  utilizing  27b 
grid  points  and  279  plate  bending  elements.  Calculated  transient  forces 
were  used,  and  a  constant  physical  mass  ot  water  ot  approximately  conical 
shape  circumscribing  the  blade,  was  taken  as  added  effective  mass.  One-hall 
of  the  imaginary  cone  of  water  was  actually  used  to  simulate  full  cavita¬ 
tion.  Transient  response,  in  the  form  oi  measured  stresses  at  a  gage  near 
the  root  of  the  blade,  was  plotted  versus  time  for  the  period  during  which 
the  blade  enters  the  water  and  reaches  maximum  submergence. 

The  SES-IOOB  propeller  was  selected  for  this  transient  response  analy¬ 
sis  because  of  t lie  availability  of  geometrical  and  test  data.  Due  to  the 
developmental  nature  of  this  application,  however,  and  to  time  and  cost 
constraints,  the  idealization  constructed  was  rather  coarse,  consisting  of 
9b  grid  points  and  nine  C1HE.Y2  isoparametric  three-dimensional  elements. 

A  detailed  study  of  this  blade  would  require  a  finer  idealization,  particu¬ 
larly  so  because  it  contained  an  annex,  or  discontinuity,  in  the  pressure 
surface  curvature.  The  exercise  of  engineering  judgment  attempted  to  com¬ 
pensate  for  this  coarseness  in  constructing  the  idealization,  in  which, 
for  instance,  thicknesses  at  grid  points  do  not  all  coincide  exactly  with 
those  on  the  actual  blade.  One  ot  the  constraints  on  placement  of  element 
boundaries  was  to  require  that  one  set  ot  boundaries  coincide  with  water 


lines  at  successive  assumed  stages  of  submergence  of  the  blade.  This 


constraint  would  not  be  necessary  with  a  fine  idealization. 

A  view  of  the  blade,  looking  in  the  forward  direction,  is  shown  in 
Figure  1.  A  series  of  water  lines  at  various  stages  of  submergence  was 
constructed  from  data  reported  by  Lewis*’  and  drawn  on  Figure  1.  An  ideali¬ 
zation  was  constructed  with  two  of  the  element  boundaries  approximately 
coinciding  with  two  water  lines,  and  this  is  shown  in  Figures  2  and  3.  In 
these  figures  the  water  lines  are  shown  slightly  curved  because  the  view 
of  the  propeller  is  flattened.  The  location  of  the  strain  gage  from  which 
Allison's  transient  response  data"1  were  obtained  is  also  shown  in  Figure  2. 
For  the  first  arbitrarily-defined  stage  of  submergence,  elements  7  through 
9  are  submerged;  for  the  second  stage,  elements  4  through  9  are  submerged; 
and  for  the  third  stage,  all  elements  are  submerged.  The  fluid  mass 
idealized  for  this  case  was  of  approximately  the  same  size  and  shape  rela¬ 
tive  to  the  blade  as  for  the  P3604  vibration  analysis,  but  the  fluid  for 
the  SES-100B  was  divided  into  only  two  layers.  An  outline  of  one  face  of 
the  fluid  mass,  showing  the  fluid-structure  interface,  is  shown  in  Figure  4. 
The  entire  fluid  mass  is  shown  in  Figure  5 • 

The  operating  conditions  for  which  this  calculation  was  carried  out 
were  the  design  conditions  of  83  knots,  1870  rpm,  and  33-degree  blade 
pitch.  Blade  pressures  are  discussed  by  Allison’  and  Lewis, h  both  of  whom 
considered  the  pressure  at  each  point  to  be  applied  in  three  stages.  First 
was  an  impact  pressure,  computed  according  to  a  theory  by  Pierson,7  fol¬ 
lowed  by  an  entry  pressure  predicted  from  a  theory  by  Wang.*’  Finally,  a 
steady  value  was  reached,  the  value  predicted  from  a  theory  by  Vim.  The 
same  procedure  was  followed  here  to  calculate  forces  to  apply  to  grid 


points  73,  75,  47,  49,  21,  and  23.  These  are  mid-side  nodes  and  are  so 
situated  that  force  is  applied  to  each  when  the  propeller  reaches  the  mid¬ 
point  of  its  stage  of  submergence.  The  pressure  versus  time  plot  at  each 
of  these  grid  points  was  derived  from  data  given  by  Lewis. ^  These  plots 
are  presented  in  Figure  6  and  show  that  these  are  essentially  pulse  pres¬ 
sures  with  no  oscillating  component. 

The  transient  response  analysis  consisted  ul  one  fluid  and  one  pro¬ 
peller  computer  run  for  each  stage  of  submergence,  a  total  of  six  runs. 
Stresses  are  printed  out  in  NASTRAN  at  the  grid  points  and  at  the  element 
centroids.  Unfortunately,  it  was  not  feasible  to  place  a  grid  point 
exactly  where  experimental  data  were  available.  Therefore,  stresses 
computed  at  three  grid  points  (22,  3b,  and  39)  which  surround  this  point 
(shown  as  Gage  1  in  Figure  2)  were  plotted  and  the  plots  are  shown  in 
Figures  7,  8,  and  9.  Static  stresses  were  also  calculated  and  plotted  on 
these  same  figures.  The  three  portions  of  the  dynamic  forces  were  averaged 
and  applied  at  the  same  points  to  obtain  these  static  stresses.  The  dyna¬ 
mic  pressures  were  also  applied,  but  with  the  steady  state  value  truncated, 
to  the  propeller  in  air,  and  the  results  are  given  in  Figures  10,  11,  and 
12.  The  plot  of  dynamic  stresses,  reproduced  from  Allison's  paper, J  is 
presented  here  as  Figure  13.  Natural  vibration  frequencies  for  the 
SES-100B  propeller  at  the  three  stages  of  submergence  were  also  calculated 
and  are  presented  in  Table  2. 

DISCUSSION  OF  SES-100B  ANALYSIS  RESULTS 

The  primary  purpose  of  this  analysis  was  to  verify  the  computational 
procedures  developed  and  to  produce  results  that  appear  reasonable  and  that 
provide  some  insight  into  the  behavior  ol  such  propellers.  It  was  also 
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desired  that,  despite  the  coarseness  of  the  idealization,  some  meaningful 


comparisons  could  be  made  with  experimental  data. 

The  plots  of  transient  stresses  in  the  propeller  in  water  (Figures  7, 
8,  and  9)  show  that  the  peak  stresses  at  the  three  locations  plotted  range 
from  21501.)  to  32000  psi.  The  measured  peak  stress  shown  in  Figure  13  is 
about  31000  psi.  This  apparent  close  agreement  is  probably  fortuitous. 

The  stress  amplitudes  shown  in  the  analytical  plots  vary  with  a  predominant 
frequency  equal  to  the  lowest  natural  frequency  of  the  propeller,  as  veri¬ 
fied  by  the  natural  frequency  analysis.  What  appears  to  be  a  similar  pre¬ 
dominant  vibrational  component,  with  nearly  the  same  frequency,  is  dis¬ 
played  in  the  experimental  plot  (Figure  13).  Substantial  amplitude  damping 
is  evident  in  the  experimental  plot,  however,  even  within  each  complete 
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period  of  vibration.  Both  the  analytical  and  experimental  plots  indicate 
that  the  blade  vibrates  through  several  complete  cycles  before  becoming 
fully  submerged.  Thus,  pressure  due  to  submergence  in  the  later  stages 
are  applied  when  the  existing  calculated  amplitudes  of  deflection  are  in 
error,  due  to  the  lack  of  damping.  Furthermore,  the  calculated  phase  could 
be  in  error  due  to  the  probable  slight  inaccuracy  of  the  calculated  natural 
frequencies.  Calculated  response  due  to  these  later  water  entry  pressures 
may  therefore  detract  from,  rather  than  improve,  accuracy  of  the  overall 
dynamic  response.  The  calculated  stress-time  plots  show  no  apparent  reac¬ 
tion  to  these  later  stage  pressures,  however.  This  is  evidently  because 
the  last  part  of  the  propeller  to  enter  the  water  is  the  portion  close  to 
the  hub  where  lower  pressures  are  generated  at  impact  (due  to  the  smaller 

*t)amping  was  not  accounted  for  in  the  analysis,  hut  could  be,  as  discussed 
later  . 
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radii  and  hence  lower  linear  velocity  resulting  from  rotation)  and  impact 
pressures  produce  smaller  bending  moments  due  to  the  smaller  effective 
moment  arms.  The  plot  of  calculated  stress  versus  time  is  probably 
meaningful  quantitatively  only  to  the  first  peak  if  damping  is  not  ac¬ 
counted  for.  Fortunately,  this  portion  of  the  curve,  or  more  particularly 
the  value  of  this  peak,  is  of  primary  interest. 

Two  other  features  of  these  curves  are  significant.  First,  peak 
stress  values  are  practically  the  same  for  the  in-water  and  the  in-air 
cases.  Second,  these  peak  values  are  on  the  order  of  two  to  three  times 
the  static  values.  This  suggests  that,  for  this  approximately  pulse-type 
pressure  application,  the  blade  behaves  somewhat  like  a  weight  suddenly 
suspended  by  a  spring:  the  weight  oscillates  about  the  static  equilibrium 
point  with  an  amplitude  double  the  static  deflection.  However,  the  ratios 
of  the  dynamic  to  static  stresses  in  the  propeller  blade  probably  vary  at 
different  locations  because  of  the  complicated  shape.  For  this  particular 
example,  a  Fourier  series  representation  of  the  assumed  pressure  pulses 
apparently  would  not  contain  any  strong  components  near  any  of  the  natural 
frequencies  in  water  or  in  air.  The  damping  effect  of  water  no  doubt 
serves  to  decrease  peak  dynamic  stresses  slightly,  but  the  main  effect  of 
submersion  in  water,  apart  from  transmitting  force  to  the  blade,  appears 
to  be  in  altering  its  natural  frequencies.  If  a  detailed  study  of  a  pro¬ 
peller  were  to  be  made,  it  would  probably  be  wise  to  obtain  stresses  at 
all  grid  points  (at  least  all  the  points  on  one  surface)  because  a  rea- 
sonance  magnification  might  occur  locally,  and  hardly  show  up  at  the  root 
of  the  blade. 
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The  fundamental,  or  lowest,  natural  frequencies  of  the  blade  in  air 


and  in  water  may  be  derived  from  the  stress  versus  time  plots  since  the 
time  record  is  probably  long  enough.  All  three  plots  for  the  in-air  case 
display  frequencies  of  about  650  Hz,  slightly  higher  than  the  calculated 
natural  frequency  of  b24  Hz  shown  in  Table  2.  The  plots  indicate  that  in 
water  the  blade  rings  at  a  fairly  constant  frequency  of  about  51U  Hz,  which 
is  close  to  tlie  natural  frequency  of  5'4  Hz  calculated  for  the  blade  under 
its  first  stage  of  submergence.  During  must  of  the  time  covered  by  the 
plots,  however,  the  blade  is  fully  submerged,  and  one  would  expect  the  fre¬ 
quency  to  fall  to  about  42U  Hz.  A  smaller  initial  stage  of  submersion,  for 
which  the  fundamental  frequency  is  higher,  might  produce  a  corresondingiy 
higher  ringing  frequency.  In  this  case,  however,  the  initial  load  would 
also  be  smaller,  being  applied  to  the  smaller  area  of  the  blade,  and  might 
not  predominate  as  it  does  in  the  present  case. 

Damping  (either  hydrodynamic  or  structural)  could  readily  be  accounted 
for  in  this  analysis,  but  the  values  of  damping  terms  would  have  to  be  de¬ 
termined  by  the  user,  probably  empirically.  All  damping  would  be  accounted 
for  in  tlie  structural  computations.  Since  damping  terms  can  be  derived 
from  the  reJative  amplitudes  of  successive  cycles  of  a  freely  vibrating 
structure,  plots  such  as  Figure  13,  in  which  the  attenuation  of  amplitude 
after  initial  response  is  clearly  illustrated,  could  probably  be  used  tor 
this  purpose.  The  principal  effect  of  damping  on  peak  calculated  stresses 
would  probably  occur  when  two  or  more  substantial  dynamic  forces  are 
applied,  separated  by  time  intervals  ot  the  order  of  one-half  a  fundamental 
period . 
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It  would  be  pointless  to  carry  the  analysis  further  in  time  to  cover 
the  emergence  of  the  blade  from  the  water,  since  the  dynamic  forces  applied 
at  that  time  would  result  essentially  from  releasing  the  relatively  small 
steady  state  submergence  pressures.  These  forces  would  be  much  less  severe 
than  those  at  entry. 

It  has  been  assumed  in  this  analysis  that  dynamic  behavior  of  a  pro¬ 
peller  can  be  adequately  described  by  analysis  of  a  single  blade,  fixed 
along  its  base.  Errors  due  to  this  assumption  may  be  significant  in  some 
cases,  particularly  for  propellers  without  relatively  massive  hubs.  The 
errors  could  be  introduced  in  two  ways.  First,  flexibility  of  the  hub 
violates  the  assumption  of  complete  fixity,  leading  .o  lower  natural  fre¬ 
quencies  than  calculated.  This  effect  was  noted  by  Allison5  in  connection 
with  the  SES-1UUB  propeller  and  with  tests  of  a  single  blade  mounted  in 
a  test  jig.  Second,  hub  flexibility  and  lack  of  great  mass  may  permit 
a  significant  degree  of  dynamic  interaction  among  the  blades.  It  is  also 
possible  that  in  some  cases  significant  interaction  may  occur  through 
the  water  between  the  blades. 

The  computational  procedures  followed  here  could  no  doubt  be  simpli¬ 
fied  and  automated  to  some  extent  with  the  aid  of  further  executive  module 
programming. 

CONCLUSIONS 

L.  Rational  procedures  have  been  presented  for  calculating  the  natural 
vibration  modes  and  transient  response  of  partially  or  fully  submerged, 
supercavitating  or  conventional  marine  propellers.  Cost  of  these  analyses 
should  be  moderate,  and  they  require  only  use  of  a  standard  NASTRAN  program. 
2.  Both  the  vibration  modes  and  the  transient  response  analyses,  particu¬ 
larly  when  progressive  submergence  is  treated,  need  more  verification. 


3.  For  the  essentially  pulse-shaped  press ure-time  profile,  as  used  here  in 
the  SES-100B  propeller  example,  peak  root  stresses  were  predicted  equally 
well  whether  the  blade  was  considered  in  air  or  in  water.  This  may  not 
have  been  true  of  stresses  elsewhere  on  the  blade.  Also,  it  may  have  been 
by  chance  that  no  resonances  of  the  blade  were  excited  by  components  of 
the  assumed  pressures.  If  time-dependent  pressures  with  oscillatory  com¬ 
ponents  were  applied,  it  is  expected  that  peak  stresses  would  be  affected 
by  submergence  in  water. 

A.  Calculated  stresses  indicated,  as  did  published  experimental  results, 
that  blade  root  stresses  due  to  pulse  pressures  can  be  two  or  more  times  as 
high  as  calculated  stresses  based  on  the  same  pressures  applied  statically. 
b.  Detailed  studies  of  stresses  and  deflections  throughout  typical  pro¬ 
peller  blades,  utilizing  finer  finite  element  meshes  would  be  useful  in 
providing  more  information  about  response  of  typical  propellers  under 
various  conditions. 

6.  When  damping  was  not  taken  into  account,  transient  response  after  about 
the  first  half  of  a  fundamental  vibration  period  was  quite  unrealistic, 

but  peak  stress  values  were  apparently  not  affected.  Damping  constants 
need  to  be  established,  perhaps  through  an  empirical  approach.  The 
apparent  failure  of  the  blade  ringing  frequency  to  fall  promptly  to  the 
fully  submerged  value  needs  further  investigation. 

7.  It  would  be  worthwhile  to  simplify  and  automate,  whenever  possible, 
the  procedures  used  here. 
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Figure  6  -  Time-Dependent  Blade  Pressures 
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Figure  Q  -  Calculated  Radial  Stress  in  Progressively  Suhmerg  i  iu'.  Propeller  Blade 
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Figure  10  -  Calculated  Radial  Stress  at  Grid  Point  22  of 
Propeller  Blade  in  Air  with  Applied  Dynamic  Forces 
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TABLE  1  -  CALCULATED  NATURAL  FREQUENCIES  OF  P3604  PROPELLER  BIADE 


Frequency 


V ibrat ion 
Mode,  N 


Calculated 

Frequenc ies 

In  Air 

In  Water 

461 

225 

1228 

671 

1680 

774 

2397 

1171 

2613 

1315 

_  .  in  Water 

Ratio,  — r — — - 

in  Air 


TABLE  2  -  NATURAL  FREQUENCIES  OF  SES-100B  PROPELLER  BLADE 


V ibrat ion 
Mode,  N 

Calculated 

Frequenc ies 

i  n  Air 

In  Water 

Stage  1 

Stage  2 

Stage  3 

1 

624 

524 

425 

420 

2 

1399 

1037 

924 

916 

3 

1679 

1316 

1206 

1158 

, 

2488 

1893 

1  581 

5 

2612 

204  5 

1815 

1797 

APPENDIX  -  INPUT  DATA  PREPARATION 


Input  required  for  the  vibration  and  transient  response  analyses 
in  this  report  wifi  be  described  using  the  SES-1U0B  calculations  as 
examples.  Inputs  described  will  be  those  required  for  the  fluid  II  and 
L  matrix  runs,  propeller  vibration  runs  in  water,  and  propeller  transient 
response  runs,  with  progressive  submer<'ence  of  the  blade.  The  fluid  matrix 
is  the  same  for  either  the  vibration  or  the  transient  response  runs. 

NASTRAN  input  decks  consist  of  control  cards,  an  Executive  Control 
deck,  a  Case  Control  deck,  and  a  Bulk  Data  deck.  Formats  for  all  cards 
used  here  are  defined  in  the  NASTRAN  User's  Manual  .  The  DMAP  ALTER  in¬ 
structions  contained  in  the  Executive  Control  Deck  and  the  checkpoint 
dictionaries  apply  only  to  NASTRAN  Level  17. 

FLUID  MATRIX  CALCULATION 

To  obtain  a  series  of  fluid  H  and  L  matrices,  one  for  each  stage 
of  submersion,  one  static  stress  computer  run  is  made  for  each  stage. 

These  are  basically  static  stress  (Rigid  Format  1)  runs,  utilizing  CIUEXJ 
elements.  For  convenience,  all  grid  points  should  be  assigned  numbers 
higher  than  those  to  be  assigned  to  any  grid  points  of  the  propeller. 

The  element  and  grid  point  data  are  assembled  in  the  usual  way,  with  the 
following  exceptions  and  notations. 

1.  All  degrees  of  freedom  ot  all  grid  points  on  the  outer  boundary, 
except  for  those  on  that  part  of  the  fluid-structure  interface  surface 
that  are  to  be  ii.  contact  with  the  propeller  during  the  current  stage,  are 
constrained  with  SPC  cards. 

2.  All  but  one  translational  degree  ot  freedom  (components  1,  ,  or 

3)  ot  ail  other  grid  points  are  constrained. 


3.  Material  properties  are  input  on  a  MAT1  card.  The  value  of 
the  shear  modulus  entered  is  the  reciprocal  of  the  mass  density  of 
water,  1/p,  for  the  present  case  in  which  water  contacts  only  one  side  of 
the  blade.  Young's  modulus  is  a  very  large  factor,  say  1U“U,  times  the 
shear  modulus.  No  value  is  input  for  Poisson's  ratio,  so  that  the  NASTRAN 
program  will  compute  a  value. 

A.  Unit  pressure  is  applied  on  the  element  faces  jf  the  fluid  - 
structure  interface  which  are  to  be  in  contact  with  the  propeller  in 
the  present  stage.  Input  is  by  means  of  PL0AD3  cards.  The  unit  pressure 
can  be  either  positive  or  negative,  as  the  sign  will  not  affect  the  final 
result . 

5.  All  grid  points  in  contact  with  the  propeller  during  this  stage 
are  listed  on  ASET1  cards. 

The  Executive  Control  and  Case  Control  decks,  and  two  additional  cata¬ 
log  instructions  which  must  be  included  in  the  control  cards,  are  shown 
in  Figure  A.l.  The  two  instructions  identify  two  permanent  files, 
RMATR1XSTEP1 ,  and  L0ADSANDTABLES1 ,  which  will  contain  the  H  matrix  and 
other  information,  respectively. 

For  the  second  and  succeeding  stages  ot  submersion,  the  only 
physical  change  is  that  a  new  set  of  grid  points  comes  into  contact  with 
tlie  propeller  blade.  This  set  will  generally  include  all  the  gridpoints 
included  in  the  previous  stage  of  submersion.  For  each  succeeding  run, 
then,  the  grid  points  newly  in  contact  with  the  propeller  must  be  added  to 
the  set  of  ASET1  points  and  removed  from  the  SPC  set  lor  which  all  degrees 
of  freedom  are  constrained.  The  element  faces  newly  in  contact  with  the 
propeller  must  be  added  to  the  set  of  loaded  element  t  rices  ,  and  the  EATAl.ut, 
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cards  must  be  changed  to  assign  new  names  to  the  H  matrix  and  other 
data.  The  names  used  here  for  the  second  step  were  HMATKIXSTEH2  and 
LUADSANDTABLES2. 

The  printout  of  the  "PC"  matrix,  or  list  of  equivalent  forces,  is 
shown  in  Figure  lb.  It  was  necessary  to  obtain  this  matrix  before  con¬ 
straints  were  applied  in  order  to  obtain  all  three  components  tor  each  grid 
point;  consequently,  the  zeros  for  all  other  unloaded  components  were 
also  printed  out.  Each  group  of  tiiree  non-zero  numbers  represents  the 
three  components  of  force  at  a  fluid  grid  point  in  contact  with  a  propeller 
grid  point.  The  numbers  are  given  in  the  order  in  which  the  fluid  grid 
points  are  numbered.  These  values  are  to  be  entered  as  1.'  matrix  terms 
on  liMlG  cards  on  which  the  matching  fluid  and  propeller  interlace  grid 
points  will  also  be  listed.  The  procedure  tor  entering  these  terms  on 
cards  is  best  understood  by  examining  the  example  to  be  presented  later. 

VIBRATION  ANALYSIS 

The  Executive  Control  and  Case  Control  decks  tor  the  vibration  analy¬ 
sis  oi  a  submerged  propeller,  and  two  ATTACH  instructions  which  must  he 
included  with  the  control  curds  ,  are  given  in  Figure  A. 2.  The  I  ile  names 
on  these  two  cards  identity  the  il  matrix  and  other  data  which  have  turn 

created  by  a  previous  run.  In  this  example,  the  names  identify  the  third, 

or  1 ui i y-s uhme rged  ,  case.  The  ALTER  routine  in  the  Executive  i out  til  Hi  .  . 

contains  two  partition  (I’AKTN)  i  ns  t  I  in  t  i  oils  ,  one  o!  whith  contains  tin 

matrix  name,  I’VLCI.  That  name  was  arbitrarily  hoseti  to  identity  a 
partitioning,  >  ul  umn  matrix  ul  l  lie  same  name  t  ' '  he  input  in  t  tie  hull-  dot.,, 
file  hulk  data  is  set  up  has  i  <  a  1  1  v  as  a  pi  ope  1  let  vihi.it  nm  ion,  will 
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t  lie  t  id  low  i  ng  a dil  i  t  1  ons  : 


1.  All  fluid  points  in  contact  with  the  propeller  (in  this  fully- 
submerged  case,  all  fluid-structure  interface  points)  are  listed  as 
scalar  points  on  SP01NT  cards  in  this  propeller  run. 

2.  Constraints  on  the  zero  component  of  these  scalar  points  are 
input  on  SPC  cards . 

3.  Two  column  matrices,  PVECD1  and  PVEC3 ,  are  input.  In  tills 
example  the  total  number  of  unconstrained  degrees  of  freedom  is  96x6, 

or  576  for  the  propeller,  plus  40  for  the  scalar  points,  or  616.  Dimen¬ 
sions  of  tlie  two  column  matrices,  therefore,  are  1x61b.  All  entries  in 
PVECUl  are  zero,  and  all  except  the  last  4(J  (assuming  the  SPOINT  numbers 
are  higher  than  all  propeller  grid  point  numbers)  of  the  PVEC3  entries  are 

zero.  The  remaining  4U  entries  are  1.0. 

T 

4.  The  L  matrix  terms  are  input  on  DM10  cards.  In  this  case, 
aUxi  terms  are  input. 

TRANSIENT  RESPONSE  ANALYSIS 

The  procedures  tor  preparing  the  additional  bulk  data  input  for  the 
vibration  analysis  are  essentially  the  same  as  those  for  the  transient 
response  case,  which  will  be  illustrated  in  more  detail. 

The  Executive  Control,  Case  Control,  and  Bulk  Data  decks,  and  the 
tour  additional  control  cards,  for  the  first  stage  of  submersion  of  the 
transient  response  case  are  shown  in  Figure  A. 4.  Of  the  four  control  cards, 
the  two  ATTACH  instructions  retrieve  the  appropriate  H  matrix  for  this  step, 
as  in  the  vibration  calculation.  The  REQUEST  and  CATALOG  instructions 
effect  the  storage  of  the  checkpoint  tile  named  SUBSTPl,  required  for  the 
restart.  The  bulk  data  input  consists  of  the  same  data  as  tor  a  straight¬ 


forward  transient  response  analysis,  and  additional  SPOINT,  SPC,  DMIG,  and 
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DM1  data  similar  to  that  for  the  vibration  in  water  case.  The  following 
comments  apply  to  the  additional  data; 

1.  All  fluid  grid  points  to  come  into  contact  with  the  propeller 
during  the  course  of  the  analysis  (not  just  those  in  contact  in  the  current 
stage  of  submersion)  must  be  listed  as  scalar  points  on  SP01NT  cards,  and 
also  on  SPC  cards  because  scalar  points  cannot  be  added  after  solution  of 
the  problem  has  begun. 

2.  PVECD1  and  PVEC1  matrices  have  the  same  dimensions  for  all  stages, 
in  tills  example  lx  bib.  The  number  of  "l.U"  terms  in  PVEC1 ,  however,  in 
this  case  18,  must  be  equal  to  the  number  of  fluid  points  in  contact  with 
the  propeller  for  the  current  stage,  and  the  placement  of  these  terms  with¬ 
in  the  "scalar  point  section"  of  the  PVEC1  column  ( t he  last  40  numbers) 
must  correspond  to  the  positions  of  the  18  contact  points  within  the 
"scalar  point  section."  All  degrees  of  freedom  of  propeller  grid  points 
and  scalar  points  are  represented  in  the  PVEC1  column  in  numerical  order. 

In  this  example  the  fluid  points  were  numbered  so  that  the  portions  ol  the 
blade  with  the  highest  numbered  points  contact  t lie  water  first;  thus  any 

"l.U"  entries  in  PVEC1,  or  PVEC2 ,  etc.,  are  the  last  terms  in  the  column. 

T 

3.  The  L  matrix  input  by  DMIO  cards  tor  this  stage  was  made  up  from 
the  corresponding  "Pt;"  matrix  printout  from  the  first  stage  1 luid  run, 
given  in  Figure  A. 2.  The  current  run  is  terminated  when  the  time  covered 
by  the  TSTKP  card,  that  is,  number  of  steps  times  length  of  each  step  (in 
this  case  0.9093x10  ^  seconds),  lias  elapsed. 

To  continue  the  analysis  for  the  second  stage  of  submersion,  a  restart 
run  is  carried  out,  with  the  following  changes  to  the  deck: 


1.  The  two  ATTACH  instructions  and  the  CATALOG  instruction  described 
tor  the  first  run  are  replaced  with  the  following  control  cards: 

ATTACH, 0PTP,SUBSTP1,ID=CSRV. 

ATTACH, INPT,HMATR1XSTEP2 , 1D=CSRV  . 

ATTACH, UT1,  LOADSANDTABLES2 , ID=CSRV  . 

CATALOG , NPTP , S  UBSTP2 , 1 D=CSRV . 

2.  The  checkpoint  dictionary  produced  by  the  previous  run  is  added  to 
tlie  Executive  Control  deck.  Mso,  PVEC1,  which  appears  in  one  of  the 
PARTN  instructions,  is  changed  to  PVEC2.  The  following  lines  are  added 
to  the  previous  ALTER  routine  and  left  for  the  remaining  restart  runs: 

ALTER  1 4 J  ,  1 4 i 

CHKPNT  K2PP,  K2DD,M2DD,  b2DL), M1U>, KDI)  S 
ALTER  149,  1  9 

I.  The  TSTEP  card ,  TslEP=2i,  in  the  Case  Control  deck,  is  replaced  with 
I'STE  P  =  2  2  . 

*.  I  he  Hulk  I  Li  t  a  deck  consists  ot  the  tol  lowing  cards: 

i.  A  delete  (,  )  card  which  eliminates  the  previous  DM It;  cards. 

T 

h.  c.irils  ,  which  represent  the  I.  matrix  lor  stage  two. 

.  .  i  >M  I  cards  tor  a  new  column  matrix,  PVEC2.  Dimensions  are  t  lie 
same  as  tor  I’VECl,  but  the  number  ot  "1.0"  terms  is  changed  to 
represent  the  new  set  ot  points  in  contact  with  the  propeller. 

In  this  case,  the  last  2b  terms  are  now  1.0. 
d.  A  new  IS  Ill  1 '  Hulk  Data  card  ,  identified  as  TSTEP  number  22, 
to  correspond  to  the  Case  Control  TSTEP  card.  This  card 
should  contain  the  correct  number  ot  time  steps  so  that  the 
total  time  defined  will  just  cover  the  current  stage  ot 
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submersion.  In  this  second  stage,  28  steps  were 
used . 

To  continue  with  the  third  and  any  succeeding  steps,  the  procedures 
lor  the  second  stage  are  followed.  The  REQUEST  and  CATALOG  control  cards 
and  the  CHKPNT  Executive  Control  deck  card  are  not  required  in  the  final 
run . 


N  A  S  T  R  A  N  executive  control  deck 


10  SES  PROPELLER  F  L U 1 0  IDEALIZATION 
APP  DISPLACEMENT 
SOL  1,0 

01  AG  1,0,14,22 
TIME  60 
ALTER  97 

MATPRN  KAA, ,,,//$ 

OUTPUT  1  ,,,,//C,N,-l  S 
OUTPUT  1  KAA,,,,//  $ 

JUMP  LBL7  S 

ALTER  112 

MATPRN  PG,ASET,,,//  5 

OUTPUT2  PG , EQE  X  t  N,  ASET,,//C,N,0/C,N,11S 

TABPT  EQEXIN, ,  ,,//  8 

JUMP  MEL  8 

ALTER  141 

LABEL  MEL  s 

EQUIV  PG, PGG  8 

ALTER  143 

JUMP  L8L0E°  % 

ALTER  165 

JUMP  FINIS  < 

ENDALTE0 

CENO 


CASE  CONTROL  OEC< 


TITLP  =  JfcS  PROPELLER  PLUIO 
SUBTITLE  =  ADDED  MASS  T  EPMS 
SPC  =  11 
LOAD  =  1 

DISPLACEMENT  -  all 
maxl infs  -  5 000 00 

BEGIN  BULK 

Ei^nri'  A.l  -  Fluid  Run  Partial  Input 
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SES  SPELLER  FLJID 
4  TOE  0  -ft S$  TEPHS 


TME  OFnSITy  OF  THIS  T  » I  <  IS  * .  *.  7  PERCENT 


ATTACH. 1NPT.HM»TRIXSTEP3«!0«CSRV, 
ATTACH. UT1 .L0ADSANDTABLFS3. ID-CSRV. 


**  A  ->  T  *  A  M  feXFCUTlWt 


C  0  M  T  O  0 


DECK 


ID  V  I  R  R  A  T  ION  jti  PROP  It,  «4T^ 

APP  DISPLACEMENT 
SOL  3  > C 

DIAG  i,  1,14,22 
TIME  61 
ALTER  53 

OPO  or  NAM  ICS , jP  L , SIL , USE T/ GPlD , S ILO , USE  TO ,,,,,,, EE  0, E QOYN/ V ,N , 

LUbET/V,N,LUSETG/V,N,NOTFL/V,N,NODLT/V,N,NOP$UL/VfN,NOFRL/V, 

N,NONlFT/V,N,NOTRL/V,N,NO££D/C,N,  /V , N, NOUE  2 
save  NOF.f  O.LUSETO  $ 

MTRXIN  CAS£CC,MA7PC0L,EQCYN, , /K2QPP , M2  OP P , 02 PP/ V , N , L USET 0/ V 
N,NOK23PP/V, N,N0M2UPP/V,N, N032PP  l 
PAPTN  K2DPP, ,PV£C3/UNUSF0, ATMAT, ,/C,N,l  $ 
iNPUir 1  /, , , ,/C,N,-l  £ 

INPUTII  /MMAT,,,,/  i 

SOLVE  HMAT , AIMA i /HI AT/C,N, 1  £ 

MPYAQ  ATMAT, HI  AT , MGG/NGGTEM/C , N, 1  2 

PARTN  MGGTEM,PVECC1,/MGGNEH,  ,  ,/C,N,-l/C,N,2/C,N,6  $ 

PARAM  //C,N,MPY/V,N,PDUM/C,N, l/C,N,-l  s 
EQUlV  MGGNEW,MGG/ POUK  $ 

JO  MP  u I C  K  £ 

MATPRN  MGGN£h,mGG, ,  ,//  £ 

LABEL  DICK  4, 

ALTER  RE , 07 
ENOAlTER 
CHKPNI  YtS 
CENT 


lASE  CONTROL  DECK 


TITLE=SES  PROP  VIBRATION 

SUBTITLE  -  PARTIALLY  SUBMERGED  PROPELLER 

LABFL  *  INVERSE  POWER  METHOO 

ECHO  =  BOTH 

K2PP  =  STEF 

me  THon=m 

SPC=1C 

DISPLACEMENT  =  All 
MAXLINES  -  SIOTcJ 
BEGIN  BUL  K 


Figure  A.  i 


Propeller  Vibration  Run  Partial  Input 


Figure  A. 4  -  Propeller  Transient  Response  Run  Input  for  Stage  1 


REQUEST tNPTDt *PF. 

ATTACH tTNpTtHMATRIXSTEPI «ID»CSRV, 
ATTACM,UTl,LOAnSANOTABLPSl.!D«CSRV. 
CATALOGtNPTPtSUBSTPl .Tn»CSRV. 


NASTRAN  executive  control  deck  echo 


TO  TRANSIENT  RESPONSE  OF  SES  PROPELLER  IN  WATER 
APP  DISPLACEMENT 
SOL  9,0 

OtBG  1,8,14,22 
TTME  60 
ALTER  72 

OPO  DVNAMTCS.ppl ,STL,USET/GPLO,5ILO,USETO,TFPUOL,nLT, ,,NLFT,TRL , , 
E(50YN7V,N,LUSET'/v,N,LUStTl>/ V, N ,N0 TEL / V , K , NO 0L  TV  V, N , NOPSOL / V , 
N.NOFRL/V ,N,N9NLFT/V,N,NOTRL/ V , N, NO  EE D/ C , N , / V , N , NOUE  f 
SAVE  LUSETB,NOOLT,NONLFT ,notrl,noue  S 

MTRYIN  P  AS  EEC  »H  ATPOOL  ,  EQDTN  ,  ,TPpOOL/K20P  p,M2DPP,B2PP/V,N,LI!SFTP/V,N, 
N0K20PP/V ,N,NOH20PP/V,N,NOe2PP  t 
P  A  VF  N(lK2r'PP,NOM2HPp,NOp2RP  X 
rARTM  K20PP, ,PVEC1/UNUSE0, ATNAT, ,/C,N,l  T 
INTUTT1  7,,  ,  ,/C,N  ,-l  i 
TNPUTT1  /HMAT, ,, ,/  $ 

SOLVE  HMAT,  AiTM»TTHI«T/C,N,  1  * 

MONAD  ATN  AT, VIAT, WGG/MGGTEMTC, N,  1  * 

EARTH  MSGt£M,PVcChi,/PG6NFM,, ,/C,N,-l/C,N ,2/C,N,f  K 
EGUJTV  liGG  new, MG'"  * 

JUMP  OTCK  * 

MATPRN  MGGNEW.MGF, ,,//  t 

LABEL  OICK  S 

ALTER  114,115 

ALTER  12*  ,129 

enoalter 

chk°nt  y  ps 

cend 
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Figure  A. 4  (Continued) 
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